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Abstract 

The new method of nonperturbative calculation of the beta-function in the 
lattice gauge theory is proposed. The method is based on the finite size 

scaling hypothesis. 

Ever since the pioneering work by Creutz [1] the approach to asymptotic 
scaling, and thus the continuum limit, was one of the central issues in studies 
of gauge theories on the lattice. Although the first results were promising, 
the lack of asymptotic scaling of physical observables has been observed in 
SU(N) gauge theories. One of the main source of the nonperturbative results 
in the gauge theories today is the Monte-Carlo (MC) lattice calculations. For 
the SU(N) pure gauge theories on lattices of size N T x N ( l MC results are the 
dimensionless functions of the bare coupling constant g (another form for the 
coupling, (3 = 2 N/g 2 , is often used). The transformation of these functions 
to physical quantities are done by multiplying them on lattice spacing a in 
the corresponding powers. The length scale L and the temperature T are 
given as 

L = AU, T={N t cl)~ 1 . (1) 

To define the physical quantities one needs a connection between lattice 
spacing a and bare coupling constant g. Such a connection is formulated in 
terms of the beta-function j3j{g) through the equation 

U(9) = ~af a ( 2 ) 

The perturbation theory gives the asymptotic expansion of the beta-function 

Pf F = -bog 3 - big 5 + 0(g 7 ), 
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where N = 2 in the SU(2) case. The differential equation (2) with (3j (g) in 
(3) leads to 
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where A^ F is the renormalization group invariant parameter (integration con¬ 
stant of Eq. (2)). Eq. (4) is known as asymptotic freedom (AF) relation. 
Using (1) and (4) one can calculate 
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The values of TJ A^ F at different N T are presented in Table 1. One observes 
rather strong dependence of T c / A^ F on N T . This means that the perturbative 
AF relation (4) does not work even on the largest available lattices. This fact 
is known as absence of the asymptotic scaling. 

It has been proposed in Ref. [3] that a deviation from the asymptotic 
scaling can be described by a universal non-perturbative (NP) beta-function, 
i.e. (3f(g) is the same one for all lattice observables and it does not depend 
on the lattice size if N a and N T are not too small. 

The following ansatz was suggested [3]: 


aAf = A (g 2 )R(g 2 ), 


( 6 ) 

where R{g 2 ) is given by (4) and A (g 2 ) is thought to describe a deviation from 
perturbative behaviour. The equation (4) has been expected at g —» 0 so 
that an additional constraint, A(0) = 1, has been assumed. The values of 
T c j A^ p can be calculated then as 


NP 
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A simple formula for the function A (g 2 ) was suggested [3]: 


(7) 


A (g 2 ) = exp 
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Parameter C 3 in ( 8 ) and a new one, T*/A^= const, have been considered as 
free parameters and determined from fitting the MC values of TJ A^ p (7) at 
different N T to the constant value T*/ A^ p . This procedure gives 


T c T/A^ p = 21, 45(14), c 3 = 5, 529(63) • 10 


-4 


(9) 
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In comparison to T c / A^ F the much weaker N T dependent values of T c / A^ p 
have been obtained, which become now close to the constant value T* /A^ p 
(9). 

Despite of the phenomenological success of the above procedure of [3] the 
crucial question, regarding the existence of the universal NP beta-function 
which does not depend on the lattice size, is not solved and remains just a 
postulate. A principal difference of our approach is that we do not assume 
the existence of the universal beta-function and take into account the finite 
size effects of the lattice. 

Usually finite size scaling (FSS) in the vicinity of a finite-temperature 
phase transition is discussed for lattice SU(N) gauge models without trying to 
make contact with the continuum limit, i.e. the scaling properties are studied 
on lattices N T x N% with fixed N T and varying N a , and the model is viewed 
as a three-dimensional spin system. On the other hand, in the continuum 
limit the FSS properties of these non-abelian models should, of course, be 
discussed in terms of the physical volume V = L 3 and the temperature T. 
On a N t x Nl lattice L and T are given in units of the lattice spacing a, 
therefore it is advantageous to introduce the dimensionless combination 


LT = (10) 

The scaling behaviour of the continuum theory emerges from the lattice free 
energy on arbitrary lattices, i.e. when varying N T and N a . 

Following [2] let us discuss briefly the FSS procedure. The singular part of 
the free energy density is described by the universal finite-size scaling function 
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where /3,y, zz are the critical indexes of the theory, the scaling function Qf 
depends on the reduced temperature t = (T — T c )/T c and the external held 
strength h through thermal and magnetic scaling fields 


gt = c t t (l + btt) 

(12) 

g h = c h h (1 + bht) 

(13) 


with non-universal coefficients q, c/ t , b t . bh still carrying a possible N T depen¬ 
dence. 
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The order parameter and the susceptibility are now obtained as derivatives 

of/ 
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Taking the fourth derivative of / at h = 0 it is easy to see that the quantity 


94 = 

is directly a scaling function 
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On a hnite lattice ^4 has the form 

/T 4 \ 

54 = 7l4 " 3 ’ (18) 

i.e. it is the normalized fourth cumulant of the Polyakov loop. 

Our approach is based on the two points: i) more conventional statistical 
mechanical definition of the beta-function and ii) FSS and the phenomeno¬ 
logical renormalization. Let us make the infinitesimal transformation of the 
lattice spacing a —> a' = ba = (1 + A b)a. Then 
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— a — = — inn 
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We get the definition of the beta-function for the lattice system 

/3 f (g) =-limh-. 

b->i db 


(19) 


( 20 ) 


Let us first consider the case when N T is kept a fixed one. Then N T can 
be absorbed in the non-universal constants in g t and gi, and we deal with 
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the usual form at the FSS as in the standard spin theory (see, for example 
[4]). The existence of the scaling function Q allows to develop a procedure to 
renormalize the coupling constant g~ 2 by using two different lattice sizes N a 
and N' a . Let us Lx the spatial size L = N a a and make a scale transformation 

a —> a = ba 

N a - K = N a /b. (21) 

Then the phenomenological renormalization is dehned by the following equa¬ 
tion 

Q(g- 2 ,N a ) = Q({g')-\N <T /b). (22) 

It expresses that the scaling function Q remains to be unchanged if the lattice 
size is rescaled by a factor b and the inverse coupling g~ 2 is shifted to ( g')~ 2 
simultaneously. Taking the derivative with respect to the scale parameter b 
of the both sides of ( 22 ) and using ( 20 ) it is easy to obtain the expression 

dg 2 = dlnQ/dlnN a 
da d In Q/dg~ 2 

The approximation of the derivative with respect to N a by the Lnite difference 
yields the formula for beta-function 
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Further we consider the formula (24) for the fourth cumulant g^g^ 2 , N a ), 
which is the scaling function directly. Fig. 1 presents the MC data for g 4 
on the lattices N T = 4; N a = 12,18,26,36 [5]. One can easy to see from 
(24) that beta-function has a zero at the fixed point 4 / g 2 = 2,299 of the 
renormalization transformation ( 22 ) in full accordance with a second-order 
nature of the deconfinement phase transition in SU(2) lattice gauge theory. 

The results of the calculation of beta-function according to formula (24) 
are presented in Fig. 2 for three sets N a = 12, N' a = 18; N a = 18, N' a = 26; 
N a = 26, N' a = 36. Although N a and N' a in the different pairs are not too 
close, one can see surprisingly the coincidence of the curves at 4 /g 2 > 4/ g 2 . 
This observation gives the hope that beta-function does not depend on the 
spatial size of lattice in the deconfinement phase. 

Next we consider fixed y = N a /N T: by varying N a , and therefore N T 
accordingly as is needed to reach the continuum limit. Rescaling N a and N T 
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by a factor b leads to a phenomenological renormalization by the following 
identity for a scaling function Q 



where g t (g 2 , N a ) is determined by (12). If we ignore the possible N T depen¬ 
dence of the coefficients c t and b tl then it follows from (25) 


t(g-\N T ) = t((g')-\N T /b). (26) 

In general the reduced temperature t = (T — T c )/T c is a complicated function 
of the coupling (3 = 2N/g 2 , which in the vicinity of the critical temperature 
T c can be approximated by [2] 


t=(P~ Pc) 
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(27) 


This approximation reproduces the correct reduced temperature in the con¬ 
tinuum limit, which is easy verified by using (4). Taking the derivatives with 
respect to the scale parameter b of the both sides of (26) and using (20) and 
(27) it is easy to obtain the expression for the beta-function 


< 3 /( 5 ) = -Bo(JVr)fl 3 --Bi(Afr)<7 £ 


where 


Bq(N t ) = 
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Then the equation (2) leads to 


aAp = exp 
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(30) 


Using (1) one can obtain the critical temperature T c . The problem only 
remains to calculate the derivative d/3 c /d\nN r in expression (29). The calcu¬ 
lation has been made for the SU(2) gauge theory by fitting the MC data for 
critical couplings 0^ G (N T ) with a spline interpolation and numerical differ¬ 
entiation of this curve. The result of the calculation is presented in Table 1. 
In comparison to T c / A^ F the much weaker dependence on N T of the critical 
temperature T c /Ap is observed. 
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N r 

A = 4/ 9c 

TJ Af 

df3 c /dN T 

TJ A L 

2 

1.880 

29.7 



3 

2.177 

41.4 

0.158 

25.22 

4 

2.299 

42.1 

0.086 

25.46 

5 

2.373 

40.6 

0.063 

25.38 

6 

2.427 

38.7 

0.045 

24.13 

8 

2.512 

36.0 

0.040 

24.24 

16 

2.739 

32.0 

0.017 

— 


Table 1: MC data for /3 C are taken from [2], The values of T c / A^ F are calculated from (4). 
Our results for T c j are obtained from (30). 
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